FARGO3D: A new GPU-oriented MHD code

We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on the physics of protoplanetary disks and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite-difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on either graphical processing units (GPUs) or central processing units (CPUs), achieving large speed-up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for CPUs, and have them translated automatically for GPUs.Fil: Benítez Llambay, Pablo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; ArgentinaFil: Masset, Frédéric S.. Universidad Nacional Autónoma de México; Méxic

Long-term and large-scale hydrodynamical simulations of migrating planets

We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems related with migration of multi-planet systems in gaseous disks, and to improve current results of migration of massive planets by including global viscous evolution as well as detailed studies of the co-orbital region during migration. We perform different tests using the public code FARGO3D to validate this method and compare its results with those obtained using a classical fixed grid.Comment: Accepted for publication in ApJ. For a movie describing the method, see https://youtu.be/66o0Z2lX8N

On type-I migration near opacity transitions. A generalized Lindblad torque formula for planetary population synthesis

We give an expression for the Lindblad torque acting on a low-mass planet embedded in a protoplanetary disk that is valid even at locations where the surface density or temperature profile cannot be approximated by a power law, such as an opacity transition. At such locations, the Lindblad torque is known to suffer strong deviation from its standard value, with potentially important implications for type I migration, but the full treatment of the tidal interaction is cumbersome and not well suited to models of planetary population synthesis. The expression that we propose retains the simplicity of the standard Lindblad torque formula and gives results that accurately reproduce those of numerical simulations, even at locations where the disk temperature undergoes abrupt changes. Our study is conducted by means of customized numerical simulations in the low-mass regime, in locally isothermal disks, and compared to linear torque estimates obtained by summing fully analytic torque estimates at each Lindblad resonance. The functional dependence of our modified Lindblad torque expression is suggested by an estimate of the shift of the Lindblad resonances that mostly contribute to the torque, in a disk with sharp gradients of temperature or surface density, while the numerical coefficients of the new terms are adjusted to seek agreement with numerics. As side results, we find that the vortensity related corotation torque undergoes a boost at an opacity transition that can counteract migration, and we find evidence from numerical simulations that the linear corotation torque has a non-negligible dependency upon the temperature gradient, in a locally isothermal disk.Comment: Appeared in special issue of "Celestial Mechanics and Dynamical Astronomy" on Extrasolar Planetary System

Dynamical friction with radiative feedback – II. High-resolution study of the subsonic regime

Recent work has suggested that the net gravitational force acting on a massive and luminous perturber travelling through a gaseous and opaque medium can have same direction as the perturber’s motion (an effect sometimes called negative dynamical friction). Analytic results were obtained using a linear analysis and were later confirmed by means of non-linear numerical simulations which did not resolve the flow within the Bondi sphere of the perturber, hence effectively restricted to weakly perturbed regions of the flow. Here we present high-resolution simulations, using either 3D Cartesian or 2D cylindrical meshes that resolve the flow within the Bondi sphere. We perform a systematic study of the force as a function of the perturber’s mass and luminosity, in the subsonic regime. We find that perturbers with mass M smaller than a few Mc ∼ χcs/G are subjected to a thermal force with a magnitude in good agreement with linear theory (χ being the thermal diffusivity of the medium, cs the adiabatic sound speed, and G the gravitational constant), while for larger masses, the thermal forces are only a fraction of the linear estimate that decays as M−1. Our analysis confirms the possibility of negative friction (hence a propulsion) on sufficiently luminous, low-mass embryos embedded in protoplanetary discs. Finally, we give an approximate expression of the total force at low Mach number, valid both for subcritical (M < Mc) and supercritical (M > Mc) perturbers

Gas dynamics in tidal dwarf galaxies : disc formation at z=0

Tidal dwarf galaxies (TDGs) are recycled objects that form within the collisional debris of interacting/merging galaxies. They are expected to be devoid of non-baryonic dark matter, since they can form only from dissipative material ejected from the discs of the progenitor galaxies. We investigate the gas dynamics in a sample of six bona-fide TDGs around three interacting and post-interacting systems: NGC 4694, NGC 5291, and NGC 7252 ("Atoms for Peace"). For NGC 4694 and NGC 5291 we analyse existing HI data from the Very Large Array (VLA), while for NGC 7252 we present new HI observations from the Jansky VLA together with long-slit and integral-field optical spectroscopy. For all six TDGs, the HI emission can be described by rotating disc models. These HI discs, however, have undergone less than a full rotation since the time of the interaction/merger event, raising the question of whether they are in dynamical equilibrium. Assuming that these discs are in equilibrium, the inferred dynamical masses are consistent with the observed baryonic masses, implying that TDGs are devoid of dark matter. This puts constraints on putative "dark discs" (either baryonic or non-baryonic) in the progenitor galaxies. Moreover, TDGs seem to systematically deviate from the baryonic Tully-Fisher relation. These results provide a challenging test for alternative theories like MOND.Peer reviewe